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Supercooled water is a metastable phase of liquid wa- 
ter below the melting temperature [1 . In this regime, the 
transition to the solid phase is irreversible once the pro- 
cess is activated. An interesting discussion recently devel- 
oped on the relationship between crystallization rate and 
the time scales of equilibration within the liquid phase 
Calculations using a coarse grained monatomic 
model of water, the mW model, suggested that equili- 
bration of the liquid below the temperature of homoge- 
neous nucleation Tjj « 225 K is slower than ice nucle- 
ation [3]. This observation has important consequences 
to a proposed theory of water anomalies, predicting a sec- 
ond critical point below Tjj where a liquid-liquid phase 
transition occurs [I] . Although it has attracted attention 
[SHE], this theory is not without problems. If the speed of 
ice nucleation is faster than liquid relaxation, the liquid- 
liquid transition would loose sense from a thermodynam- 
ical point of view, being the liquid phase not equilibrated 

Here, a 3 [is long molecular dynamics simulation of 
the TIP4P-Ew water model is presented to investigate 
the relaxation properties of an atomistic model in the 
supercooled region below Th- The length of this calcu- 
lation is one order of magnitude larger than the 350 ns 
used to study freezing with the mW model [3]. A box 
of 1024 molecules was simulated with GROMACS [9]. 
The Berendsen barostat [[TOj, velocity rescale thermostat 
[TT] and PME [T2] were used for pressure coupling, tem- 
perature coupling and long-range electrostatics, respec- 
tively. The simulation was run at 190 K and 1250 atm. 
These values are close to the estimated liquid-liquid crit- 
ical point for the TIP4P-Ew [13 , congruous with recent 
calculations on the similar TIP4P/2005 model [5J. 

In these conditions, freezing was not observed as shown 
by the timeseries of the potential energy E p (Fig. [lR)- 
Fluctuations are of the order of 0.5 kJ/mol per molecule 
with no systematic drift. It has been observed that once 
freezing is activated the energy drifts very quickly to low 
values of the potential energy, with large energy changes 
(e.g. roughly 5 and 2 kJ/mol per molecule for TIP4P at 
230 K QJ| and TIP4P/2005 at 242 K [TS], respectively). 

The time series of the density p and the tetrahedral 
order parameter Qt [E] are shown in Fig. [l^-C. They 
respectively correlate and anticorrelate with the poten- 
tial energy (Pearson correlation coefficient r = 0.69 and 
-0.86). The distributions of both p and Qt show an ap- 
preciable bump at one of the tails (see right panel of 
Fig. [TJ3-C) , suggesting the presence of a subpopulation. 
For the case of the tetrahedral order parameter, the sub- 
population emerges at values around 0.873 (red dashed 
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FIG. 1. Time series for the 3 ps trajectory. (A) potential 
energy; (B) density; (C) tetrahedral order parameter Qt', (D) 
parameter (calculated as in Ref. \TE\). Right panels 
show the probability distribution of the respective quantities. 
Dashed lines represent deviations from the mean fluctuations, 
with values of 0.972 and 0.873 for p and Qt, respectively. 

line and right side of Fig. [lp). This fluctuation is local- 
ized in a time window between 2.3 and 2.6 ps in corre- 
spondence to a decreasing of both the density and the 
potential energy. It is interesting to note that density 
subpopulations have been interpreted by some |18| as a 
signature of the aforementioned liquid-liquid transition. 

To check whether this fluctuation corresponded to an 
ice nucleation attempt, the Qq order parameter [21 HBIITU] 
was calculated (Fig. [Tj3) . In the time window between 
2.3-2.6 ps the value of the parameter is around 0.025, 
with no signs of ice nucleation. Moreover, no correlation 
with the energy was found (r = 10 -6 ). With a value of 
Qq for hexagonal ice expected to be one order of magni- 
tude larger [16] . no evidence for ice nucleation is found 
in the present trajectory. 

Finally, the oxygen mean-square-displacement (MSD) 
was calculated (Fig. [2]). At timescales shorter than 
one ns, water shows a subdiffusive behavior (dotted line 
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perimental validation is needed to clarify which proposed 
mechanism (if any) is closer to real water. 
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FIG. 2. Oxygen mean square displacement (MSD). The 
dashed and dotted lines represent a linear and a power- 
law (exponent equal to 0.1) regression, respectively. The 
diffusion coefficient extracted from the linear regime is of 
6.6 x 10" 9 cm 2 /s. Fhe g_msd function of GROMACS was 
used with 150 windows to improve statistics. 



in Fig. [2]). For larger times the system enters a diffu- 
sive regime, following the linear relationship MSD ps t 
(dashed line), with a maximum average displacement of 
3.47 nm after 3 /is. Taking into account that the molec- 
ular diameter is around 0.3 nm, water molecules have 
diffused for about 11.5 molecular diameters (the average 
box side length is of 3.14 nm). 

In conclusion, evidence is provided that the liquid 
phase of the TIP4P-Ew model is at equilibrium in the 
supercooled regime before ice nucleation. Our finding is 
in agreement with another ^s long simulation of super- 
cooled water with a 5-site model [18], suggesting that 
equilibration of the liquid phase below Tjj is a common 
feature of atomistic models. The mW model has shown 
to reproduce several properties of water, including den- 
sity and phase diagram [2U]. But the lack of hydrogens, 
and consequently of molecular reorientations |21j . might 
considerably speed up the time scales. We speculate that 
the differences in the relaxation kinetics between atom- 
istic models and the mW model are due to the lack of 
molecular reorientations in the latter. Clearly, further ex- 
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